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Relating chaos to deterministic diffusion of a 
molecule adsorbed on a surface 

O 

^N| 1 Abstract. Chaotic internal degrees of freedom of a molecule can act as noise 

and affect tfie diffusion of tfie molecule on a substrate. A separation of time scales 
^ ' between the fast internal dynamics and the slow motion of the centre of mass on 

the substrate makes it possible to directly link chaos to diffusion. We discuss the 
conditions under which this is possible, and show that in simple atomistic models 
with pair- wise harmonic potentials, strong chaos can arise through the geometry. 
Using molecular-dynamics simulations, we demonstrate that a realistic model of 
benzene is indeed chaotic, and that the internal chaos affects the diffusion on a 
graphite substrate. 
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1. Introduction 

Diffusion is always related to noise in some way. In the case of Brownian motion, and in 
most studies of diffusion, the noise is thermal in origin. Thermal noise is, in principle, 
generated by a deterministic many-particle system, but noise can also be generated 
by smaller deterministic systems if they possess the right dynamical properties. 

In the case of diffusion of large molecules adsorbed on a surface, the internal 
dynamics of the molecules can be such a source of noise. Experimental [30) and 
numerical |16| evidence exists demonstrating that rotational and other internal degrees 
of freedom affect the diffusion. In some cases clusters of molecules have been shown 
to diffuse more rapidly than single particles ^ [5S] . All this suggests that the internal 
degrees of freedom can play an important role. They may act as a finite heat bath and 
drive the diffusion. For this to happen, two important conditions must be met. First 
of all, the internal dynamics must be chaotic, so as to produce noise. Secondly, the 
internal degrees of freedom must be suitably coupled to the motion on the substrate, 
so that the noise can be transferred from one system to the other. A large body 
of work exists on the relationship between noise and chaos in finite systems, see for 
instance references [l5l l4l [3| [2| [271 [28] . 

In reference [5] it was shown that in a formally similar system, the bouncing- 
ball billiard, such a relation between internal noise and global diffusion can indeed 
be made. Under certain conditions, the system exhibits a separation of time scales 
between slow global coordinates and the fast "internal" dynamics of the particle. If 
the fast internal dynamics are chaotic it can be shown that they act as a source of 
noise for the slow global coordinates and lead to diffusion, the size of which can be 
estimated. If the internal dynamics are not chaotic, there is no noise, and hence no 
diffusion. 

In this paper we study the relation between time-scale separation, fast chaos, and 
diffusion of molecules on a substrate. We describe the molecules and their interaction 
with the substrate with realistic, classical, atomistic Hamiltonians and show how even 
pair-wise harmonic potentials can lead to nonlinearities and chaotic behaviour. We 
show that in benzene chaotic internal degrees of freedom do indeed lead to diffusion 
of the molecule on a graphite substrate, without any addition of thermal noise. 

In the following section [2] we first discuss time-scale separation in adsorbed 
molecules, and review the case of the bouncing-ball billiard. In section [3| we explain 
how the internal degrees of freedom of a molecule can produce noise through nonlinear 
dynamics, which conditions must be met, and which properties of the molecule may 
enhance the chaos. Finally, in section [31 we discuss a simple atomistic model for a 
benzene molecule, adsorbed on a graphite surface, and show some preliminary results, 
confirming that chaos in the internal degrees of freedom does indeed affect the diffusion 
on the substrate. 

2. Separation of time scales 

When a molecule is adsorbed weakly on a surface, i. e. it is not strongly distorted by 
the surface, the internal forces are stronger than the forces exerted by the substrate, 
and so the internal coordinates change more rapidly than the coordinates of the centre 
of mass on the substrate. The system has two separate time scales, a long one for 
the motion of the centre of mass on the substrate and a shorter one for the internal 
degrees of freedom of the molecule. More specifically, the equations of motion can be 



Relating chaos to deterministic diffusion of a molecule adsorbed on a surface 3 

split into those for the slow motion of the centre of mass, described by a (containing 
both position and momentum), and those for the fast internal coordinates b, 

a=F(a,b,i), (1) 

b=iG(a,b,t), (2) 

with F and G both well-behaved functions of order 1, and the time-scale separation 
parameter e a small number. 

Time-scale separation is very common in nature and exists on many different 
scales and levels of complexity, from the fairly simple orbits of planets and moons 
to the extremely complicated interaction between the day-to-day weather and the 
climate. In general, to simplify description of the slow subsystem, one wishes to 
eliminate the fast variables. In practice, however, the fast variables cannot simply be 
removed, but must be replaced with some effective behaviour. Several strategies exist 
for dealing with systems with multiple time scales. 

If the fast subsystem behaves (quasi-)periodically, the coupling to the slow system 
can be averaged out in the treatment of the slow system[55]. The coupling is replaced 
by an effective deterministic coupling term. 

If the fast subsystem is highly disordered and has a very large number of degrees 
of freedom, it can be approximated by an infinite heat bath |33) . A fast subsystem 
decays to equilibrium quickly, so that the coupling which affects the slow system can 
be approximated by equilibrium behaviour. Using projection-operator methods, such 
as described in van Kampen [34] or in Pick & Sauermann [lOj, it can be shown that, 
in general, coupling to a system consisting of many fast degrees of freedom can be 
approximated by stochastic forces and damping. 

2.1. Finite heat bath 

Here, we wish to consider fast subsystems which can are neither infinite nor (quasi- 
)periodic, but can be seen as finite heat baths. It has been suggested that chaotic fast 
degrees of freedom cannot be distinguished from noise [151 IH HI] • Therefore, we would 
like to replace the fast degrees of freedom by some suitably chosen noise. If the fast 
system is sufficiently well-behaved and correlation in it decays exponentially, the noise 
should be Gaussian white noise. In references |17[ [T51 \T7\ \7E\ it was shown through 
the use of projection-operator methods, that, under certain conditions, a reduction 
similar to the one for large systems is also valid if the slow system is coupled to a 
small number of fast chaotic degrees of freedom. Similar results were also obtained 
for different conditions in references [3l [2] through the use of the Zwanzig projection 
method, and in reference [9]. In these works the validity of linear response (or time- 
scale separation) is considered carefully. In all of these results, the fast chaos acts as 
a finite thermodynamic heat bath and leads to stochastic driving and damping of the 
slow system, with clear signatures [HI [2l] of the finite size of the fast system, and the 
finite amount of energy stored in it. 

When a dynamical system explores the entire phase space for almost all initial 
conditions, it is called mixing. If the fast subsystem has this property, the system's 
long-time behaviour is independent of the initial conditions. If the fast subsystem also 
has exponential decay of correlation, then the slow coordinates can be described by 
a Fokker-Planck equation for their probability density. In reference [271 123 explicit 
expressions were derived for the effective drift and diffusion in terms of averages over 
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trajectories of the fast subsystem for fixed slow coordinates. These expressions can be 
used to calculate the effective behaviour of the slow subsystem [51 [T] . 

Due to the mixing, the expressions for the diffusion and drift are uniquely defined, 
and cannot depend on the history of the slow subsystem. However, mixing is not 
strictly necessary. A somewhat relaxed but more complicated condition, that only 
the chaotic region of the fast subsystem be mixing and deformed continuously with 
the slow coordinates, suffices [6]. This property is much more common in physically 
motivated (Hamiltonian) dynamical systems with a small number of degrees of freedom 
than mixing itself. 

2.2. Time-scale separation and chaos in the bouncing-ball billiard 

In order to gain understanding of this kind of time-scale separation in physical systems, 
in reference [^ de Wijn and Kantz studied a simpler but still physical discrete 
dynamical system called the bouncing-ball billiard. It is a simple model for transport 
and is formally similar to transporting systems with internal degrees of freedom, such 
as molecules adsorbed at surfaces. 

The bouncing ball billiard is a two-dimensional extension of the one-dimensional 
bouncing-ball problem, which has been studied in great detail both experimentally 
[21 mi and theoretically [H [501121111] ■ It was first introduced by Matyas and Klages 
in reference Q2i- As shown in figure [H a point particle subject to a gravitational 
field g bounces inelastically (restitution coefficient a parallel to the surface and (3 
perpendicular to it) on a vibrating plate (frequency uj and amplitude A). The vibrating 
plate has a symmetric tent shape, with small slope E. 

The vertical motion of the point particle corresponds to the internal degree of 
freedom of the molecule, while the horizontal motion on the plate with a periodic 
shape corresponds to the motion of the centre of mass on the periodic substrate. 
Note that the "internal" degrees of freedom of the bouncing ball are not Hamiltonian. 
However, what is more relevant for the issues of effective diffusion is the structure 
of the phase space: as is the case for Hamiltonian systems with a small nunrber of 
degrees of freedom the system can exhibit both chaotic and periodic behaviour (is not 
mixing). One-dimensional Hamiltonian systems, on the other hand, are never chaotic. 

For small slope, the bouncing-ball billiard has a time-scale separation between 
the horizontal and vertical motion. This implies that if the vertical motion is chaotic, 
the horizontal motion should be diffusive, while if the vertical motion is periodic, 
the horizontal motion should be deterministic. This was confirmed using molecular- 
dynamics simulations. Due to the relative simplicity of the dynamics, the effective 
diffusion and drift can be evaluated (see equations (I8)-(23) of reference [6j). When 
compared to the numerical simulations, the results show that, in this model system, the 
time-scale separation theory for finite chaotic fast systems can be used to successfully 
predict the diffusive behaviour and approximate diffusion and drift. 

3. Sources of chaos and nonlinearity in atomistic models 

Through a separation of time scales, chaos in the internal degrees of freedom of a 
molecule can lead to diffusion on a substrate in the same way as the vertical chaos 
in the bouncing-ball billiard leads to diffusion in the horizontal profile of the driving 
plate. In order for the molecule to generate noise and for this to contribute to the 
diffusion, the internal dynamics must be chaotic even without external fields. When 
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Figure 1. The bouncing-ball billiard with all its parameters. The horizontal 
motion on the periodic profile of the driving plate, is formally similar to the motion 
of a particle on a substrate. The bouncing ball itself corresponds to the internal 
degrees of freedom of a molecule. 



a system is chaotic, it is exponentially sensitive to initial conditions. The rates at 
which infinitesimal perturbations of the initial conditions blow up for long times are 
referred to as the Lyapunov exponents. Given the finite resolution of any observation, 
a trajectory with such a sensitivity becomes unpredictable for long times. Seemingly 
random noise is thus generated deterministically from the infinite detail of the initial 
conditions. 

The chaotic hypothesis, put forward by Gallavotti and Cohen [T2j [H] states that 
many-particle systems, such as the ones typically studied in statistical physics, are 
strongly chaotic. For systems with fewer degrees of freedom, such as the molecules 
studied in this paper, this assumption cannot be made. We must first investigate how 
and when chaos appears in such systems. 

Only nonlinear dynamical systems may be chaotic [71126). In general, the nature 
of a trajectory in a nonlinear system depends on the initial conditions. It is not 
unusual for both regular and chaotic behaviour to occur in the same system, but in 
different parts of the phase space. 

3.1. Atomistic models 

In principle, molecules should be treated quantum mechanically. The study of chaos in 
quantum mechanical systems is a field of its own. For more details, see, for instance, 
reference [14j . In general, it is possible to treat a quantum mechanical dynamical 
system classically if the energy spacings are sufficiently small compared to the total 
amount of energy available to excite states. 

Here, we use realistic, classical potentials to describe molecules. Many such 
potentials exist in the literature, but here we choose to start from the much used 
expressions of the tripos 5.2 force field [5]. In this description, the bending and 
stretching vibrations around equilibrium positions contribute purely quadratic terms 
to the potential. Most atomistic descriptions are similar, in the sense that the pair-wise 
potentials for bond length and angles are nearly or completely harmonic. Naively, one 
would expect that pair-wise harmonic potentials cannot lead to nonlinear dynamics, 
but, as we will demonstrate below, the geometry of a molecule can still lead to 
nonlinear behaviour. 
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Figure 2. A simple atomistic model of a molecule with two atoms attached to a 
relatively heavy atom or group. 



3.2. dimer 

The model systems used here are all Hamiltonian and therefore conserve energy, as 
well as phase space volume and several other quantities. Conserved quantities and 
symmetries of the dynamics can be used to generate perturbations which do not grow 
or shrink exponentially (Lyapunov exponent zero), and therefore do not contribute 
to chaos. Perturbations in the centre of mass coordinates, for instance, grow at most 
linearly, not exponentially. 

Due to the symmetries and conserved quantities, there needs to be a certain 
minimum number of coupled degrees of freedom before a Hamiltonian system can 
display chaos. A one-dimensional dimer without external fields, for instance, cannot 
be chaotic, as it has four-dimensional phase space, and also four zero modes. Three 
zero modes are associated with perturbations in the position and momentum of the 
centre of mass, and time-translation invariance. Due to the conservation of phase- 
space volume, the fourth Lyapunov exponent must also be zero. 

In three dimensions, a dimer has a 12-dimensional phase space of which the six 
components related to the centre of mass decouple, and cannot produce chaos. The 
rotational symmetry generates four perturbations that cannot grow. Finally, there is 
time-translation invariance, combined with the conservation of phase-space volume. 
This leaves no degrees of freedom for the production of chaos. For the possibility of 
chaos to exist, therefore, the molecule must consist of at least three particles. 

3.3. Nonlinearity from the geometry 

The simplest three-particle system is a trimer, such as displayed in figure [3 For 
simplicity, we take one of the particles at the end of the trimer to be much heavier 
than the other two and this particle's position and orientation are kept fixed. The two 
lighter atoms have masses M and m respectively, while the bond lengths are R and 
r. The angles a and cj) finally describe the geometry of the bonds. In the equilibrium 
configuration R, r, (j), and a take the values i?o, tq, (po, and ao respectively, at the 
global energy minimum. Without loss of generality, ao can be taken to be 0. 

As mentioned above, we use a pair-wise harmonic potential for the bending and 
stretching. The typical frequencies associated with the bending are lower than those 
associated with the stretching, see for instance the constants in reference [5]. With 
the low frequencies come relatively large amplitudes for the same amount of energy, 
and so nonlinear terms in the angles or their time derivatives are larger than those in 
the bond lengths or their derivatives. For the purpose of demonstrating the existence 
of large nonlinear terms in the dynamics, it therefore suffices to keep the bond lengths 
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fixed. 

Let kff, and fc^ be the spring constants of the harmonic terms for (j) and a 
respectively. The total potential energy can be written as 



V 



1. 



*o)^ + ^fcatt^ 



(3) 



The kinetic energy expressed in these coordinates is more complicated, 
1 



T 



-MR^a 
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dt 



[R sin a + r sin(7r — cf) + a)] 



+ [ —[R cos a + r cos(7r — + a)] 



-MR^a^ + -m 



2\ -2 I „2l2 



2-- ^...^{R'+r')a' + r'c^- 

— 2(j)ar'^ + 2rRa{a — 0) cos(0 — n) 



(4) 



(5) 



which contains harmonic, but also strongly nonlinear terms. 

If the kinetic energy is expanded around the energy minimum, we see that all 
non-harmonic terms originate from the cos{(f) — n) factors. For 00 — t^: these higher 
order terms become fourth order, whereas for (f>Q ^ n, they are third order. For 
molecules with bond angles not equal to tt, the nonlinear terms are stronger and can 
be explored at lower energies. The geometry produces nonlinearity, even in molecules 
with pair-wise harmonic potentials, and the strength of the nonlinearity depends on 
the bond angles. 

4. Diffusion of benzene on graphite 

The considerations above suggest that one should look for chaos in molecules with a 
sufRcient number of atoms, and with somewhat complicated geometry. Molecules such 
as carbon dioxide, which has the right kind of geometry, but only three atoms, may 
not be sufficiently nonlinear for energies corresponding to room temperature. Larger 
molecules, for instance benzene, have more degrees of freedom, and are more likely to 
produce chaos. 



4.I. An atomistic model of benzene 

Consider the atomistic model of benzene shown in figure |31 The angles between 
bonds are 2tt/3, and the six heavy carbon atoms alone should have a sufRcient 
number of degrees of freedom to produce chaos. For the interaction between the 
atoms, we use the tripos 5.2 force field described in reference jHj. The bending and 
stretching contribute to the potential energy through quadratic terms around their 
energy minima. Additionally, there is also a somewhat weaker contribution due to 
torsion, which depends on the positions of 4 atoms. 

The hydrogen atoms in benzene are light, and the frequencies associated with 
various stretching and bending modes of the C-H bonds are high compared to 
frequencies related to interactions between heavier atoms. Even at room temperature, 
the quantum states of hydrogen bonds are far apart compared to the typical available 
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Figure 3. A schematic representation of the atomistic model for benzene. The 
hydrogen atoms and their bonds to the carbons are difficult to treat classically 
and have therefore been eliminated. 



energy. The associated modes are difficult to excite and tfiese bonds remain in tlie 
ground state. In otlier words, the hydrogen bonds do not participate in the chaotic 
dynamics. If an atomistic model for benzene is to be used to investigate the chaotic 
dynamics, the C-H bonds must first be eliminated. 

The hydrogen atoms can be averaged out using a mean-field approximation, by 
assuming that a hydrogen atom is, on average, in the same plane as the three nearest 
carbon atoms, i. e. the one to which it is bonded and its two neighbours. This greatly 
simplifies the torsion terms. Let r^ and pi denote the position and momentum of the 
i-th carbon atom, ordered in such a way that i and i + 1 are neighbours. Let (j^i be 
the angle between ri_i — Tj and Tj+i — r^, and /3i the torsion angle associated with the 
bonds between the i — 1-st, i-th, i + I-st, and i + 2-nd carbon atoms. The Hamiltonian 
can be written as 



H = 



1 

2(mc + m-n) 


6 


.i.. 


6 

i=i 


i=l 


'(J+l)(mod 6) - 


-r.ll 






rof 



^fc^[l-|-cos(2A)] 



(6) 



where fc^ and kp are the effective bending force constant and the effective torsion 
constant. The effective bending constant contains contributions from C-C-C, C-C-H, 
and H-C-C angles. The effective torsion constant consists of the four torsion terms for 
the combinations C-C-C-C, C-C-C-H, H-C-C-C, and H-C-C-H, 

fc0 = fcc-C-C + 2^c-c-H , (J) 

fc/3 = ^H-C-C-H + fcc-C-C-C + 2fcc-C-C-H ■ (8) 

The constants from reference 0, which will be used in this model for the rest of 



Relating chaos to deterministic diffusion of a molecule adsorbed on a surface 9 

this paper are tq = 1.47 A, k^ = 1400 kcal/ (molA^ ) , k^ = 158 kcal/(rad^mol), and 
fe^ = 5.696 kcal/mol. 

4-2. Numerical simulations of benzene 

We have performed molecular dynamics simulations of this model of benzene and 
numerically obtained the largest Lyapunov exponent, which quantifies the chaotic 
behaviour. The system was started from random initial conditions and integrated 
using the velocity- Verlet algorithm. The velocity- Verlet algorithm is symplectic and 
ensures that the conserved quantities in Hamiltonian systems remain conserved for 
finite time-step size [211 H^ . 

For a nonlinear dynamical system with a large number of dimensions it is 
impossible to determine numerically within reasonable time the set of initial conditions 
that have a specific energy. We therefore choose the initial conditions with zero total 
momentum and zero angular momentum, and on the intersection of this with the 
constant energy shell of the linearised system, with energy corresponding to a specific 
temperature. For all temperatures at which the classical description is valid, the 
Lyapunov exponents were found to be positive. At energy corresponding to room 
temperature, 32 set of initial conditions were investigated, and all of them produced 
chaos. 

4-3. Benzene on a graphite substrate 

We have extended our atomistic molecular dynamics simulations of the benzene 
molecule described before by adding an external potential to represent the coupling 
to the hexagonal lattice of a graphite substrate. We choose graphite, because of its 
weak interaction with adsorbed molecules. Each carbon atom is affected separately 
by the potential. The inter-atomic distance in the graphite was chosen to be equal to 
1.42 A and the corrugation was 25 meV. In the direction perpendicular to the surface, 
the potential was harmonic with spring constant 61.2 meV/A^. The theoretical 
expressions for the effective diffusion and drift can now, in principle, be evaluated, 
along the lines of the calculations in [6l [1] . Though the qualitative behaviour of the 
system is derived in this paper using very simple arguments, it should be noted that 
the quantitative evaluation of the theoretical expressions for the diffusion and friction 
is non-trivial, because the benzene molecule is far more complicated than the bouncing 
ball of reference [6] or the hydrogen atom of reference [1] . 

Initial conditions were again chosen randomly, with the energy of the linearised 
system corresponding to room temperature, and with zero total momentum. The 
model molecule was placed on a substrate in an energy minimum, and the system was 
integrated with time steps of size of 1.0 a.u=0.0242 fs. 

The motion of the centre of mass is displayed in figure |4] for eight different sets 
of initial conditions. All eight systems have a positive largest Lyapunov exponent. 
Though the initial conditions for the centre of mass are identical in all cases, the eight 
resulting trajectories are different. The chaotic internal dynamics therefore cause 
diffusion of the centre of mass on the substrate. 

Analysis of the time-dependence of the mean square displacement suggests that 
the diffusion is normal. However, the diffusion sometimes exhibits long jumps, which 
produces a dependence of the mean square displacement on time which only becomes 
linear for long times and for two trajectories displayed in the figure |3] the analysis is 
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Figure 4. The trajectories of ttie centres of mass of benzene molecules starting 
from identical initial conditions for the centre of mass, but different initial 
conditions for the internal degrees of freedom. The initial condition is {x, y) = 
(0, 0), {px,Py) = (0, 0). Distances are in a.u. (ao = 0.529177A). The motion on the 
substrate is diffusive, and, for sufficiently long trajectories, the diffusion appears 
to be normal. 



inconclusive. Note that in a system which is also subject to thermal noise, extremely 
long jumps could be more rare. 



5. Conclusion 



The internal degrees of freedom of molecules can be related to the diffusion of their 
centres of mass on a substrate. If the internal degrees of freedom are chaotic, they can 
be seen as a finite heat bath which generates noise that leads to diffusion of the centre 
of mass. This enhances the diffusion of a molecule with internal degrees of freedom 
on a substrate compared to diffusion with purely thermal origin. 

We have demonstrated in this paper that atomistic models of molecules with 
pair-wise harmonic potentials can produce chaotic behaviour. In a realistic model 
of benzene, chaos occurs at energies comparable to realistic temperatures. We have 
shown that the time-scale separation combined with internal chaos of the benzene 
molecule implies that the benzene molecule should diffuse on a sufficiently weak 
substrate, such as graphite without any external sources of noise. The preliminary 
numerical studies of benzene on graphite presented here confirm this. 

These results raise the possibility that the time-scale separation between the 
internal dynamics and the centre-of-mass motion could enhance the diffusion of heavy 
molecules in experiments. Within the frame-work of the time-scale separation theory, 
theoretical expressions exist for the effective diffusion of such molecules due to their 
internal degrees of freedom. However, evaluating these expressions is not a trivial 
exercise, as even the simple model of the benzene molecule on a graphite surface 
described in this paper is more complicated than the systems for which they have been 
evaluated in the past. Additionally, the interplay between the noise from the finite 
heat bath and the thermal noise from the environment will have to be investigated 
further. 
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